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Abstract 

We derive cosmological constraints on the annihilation cross section of dark 
matter with velocity-dependent structure, motivated by annihilating dark matter 
models through Sommerfeld or Breit-Wigner enhancement mechanisms. In models 
with annihilation cross section increasing with decreasing dark matter velocity, big- 
bang nucleosynthesis and cosmic microwave background give stringent constraints. 



1 Introduction 



In a weakly-interacting massive particle dark matter (WIMP DM) scenario, a DM particle 
with mass of (9(100) GeV - 0{1) TeV should have an (thermally-averaged) annihilation 
cross section of (av) ~ 3 x 10~^^ cm^/s in order to reproduce the observed DM abundance 
due to the thermal production. On the other hand, recently reported excesses of cosmic- 
ray positron [1] and electron fluxes [2l|3lll] may be interpreted as signatures of annihilating 
dark matter with fairly large annihilation cross section of order of 10~^^ - 10~^^ cm^/s 
depending on DM mass m, which is typically three orders of magnitude larger than 
the canonical value quoted above, although constraints from other observations, such as 
gamma-rays [SI El E] and neutrinos [SI El [TOl [H] are also stringent and might have already 
excluded some parameter regions. 

One way to achieve the "boost factor" of (9(10^) is to make the DM annihilation cross 
section velocity-dependent. In this case the annihilation cross section in the early Uni- 
verse is not same as that in the Galaxy or elsewhere, simply because typical velocity of the 
DM particle varies from place to place. Hence it is in principle possible that the DM has 
canonical annihilation cross section at the freezeout epoch in the early Universe reproduc- 
ing the DM abundance observed by Wilkinson Microwave Anisotropy Probe (WMAP), 
while explaining the cosmic-ray positron/electron excesses. A common mechanism would 
be Sommerfeld enhancement of annihilation cross section [121 [131 113] • If a DM interacts 
with a light particle through which it annihilates, non-perturbative effects enhance the 
annihilation cross section. The cross section is enhanced by inverse of the DM velocity, 

or f in this class of models. In the Breit-Wigner enhancement scenario on the other 
hand, DM annihilates through S-channel resonance where a particle in the intermediate 
state has a mass close to two times DM mass [151 llHl [IZ]- In this case the DM cross 
section can scale as at an earlier time or at a later time. 

In these models the annihilation cross section increases as the temperature decreases 
in the early Universe, and hence DM continues to inject high energy particles through the 
cosmic history. Therefore, it is quite non-trivial whether these models satisfy constraints 
from big-bang nucleosynthesis (BBN) and cosmic microwave background (CMB). In the 
case of velocity-independent annihilation cross section, bounds from BBN [T8l [T9l l20l 
[211 [22] and CMB [251 EH [25] were derived in previous works. In this paper, we extend 
the analysis to the velocity-dependent annihilation cross section and derive general upper 
bound on the annihilation cross section. 

This paper is organized as follows. In Sec. [2] a simple prescription for treating the 
velocity-dependence of DM annihilation cross section is described. In Sec. [3] we present 
constraints from BBN and CMB and give implications on DM models. Sec. [H is devoted 
to conclusions and discussion. 
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2 Dark matter with velocity-dependent cross section 



2.1 Models of velocity-dependent annihilation cross section 

Below we briefly give examples of DM with velocity-dependent annihilation cross section. 
After that we will explain our unified treatment for describing the cosmological effects 
from DM annihilation with velocity-dependent annihilation cross section. 



2.1.1 Sommerfeld enhancement 

A DM particle x is assumed to have an interaction with 0, which may be a scalar or 
gauge boson with coupling constant a^, whose mass is much lighter than the DM mass: 
m<^ -C m. Let us consider the DM annihilation process mediated by exchanges. If the 
mass of is sufficiently small, the 0-mediated interaction can be regarded as a long-range 
force and such an annihilation cross section receives an enhancement S compared with 
tree-level perturbative expression [26] . 

where v is the initial DM velocity in the center of mass frame. Thus the DM annihilation 
cross section is proportional to 1/f for v ^ a^. This Xjv enhancement saturates at 
V ~ m^/m. There is another interesting effect caused by the bound state formation, 
which resonantly enhances the DM annihilation rate for some specific DM mass [111 [Tl] . 
It is known that the enhancement is proportional to near the zero-energy resonance, 
and this f ~^ behavior also saturates due to the finite width of the bound state. 



2.1.2 Breit-Wigner enhancement 

In the Breit-Wigner enhancement scenario, DM particles annihilate through ^-channel 
particle exchange (0), where the mass of 0, m^, is close to 2m. The square amplitude of 
this S'-channel process is proportional to 

l-Mpoc (2) 

[y^ + by + Y 

where b and 7 are defined as = 4m^(l — S) and 7 = V^jm^ (F,^ is the decay width 
of 0), respectively. If b and 7 are much smaller than unity, we have (af) oc in the 
limit ^ max[(5, 7]. At smaller velocity, it becomes proportional to f^^. Finally, in 
sufficiently small v the cross section saturates at a constant value. 



2.2 Energy injection from DM annihilation 

Some DM models have velocity-dependent annihilation cross section as described above. 
In order to treat the effects of velocity-dependence, we phenomenologically parametrize 
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the annihilation cross section as 



{(^v) = (3) 
e + [v/vo) 

where {ctv)o is a constant, v is the (thermal-averaged) velocity of DM particle, and Vq is the 
velocity at the freezeout of DM annihilation in the velocity-independent case. Typically, 
the freezeout temperature is given by Tfo ~ m/25 ^7\, which gives Vq ~ 73/5 ~ 0.3 ; m 
our numerical calculations, we take vq = y/S/b independently of n although the freezeout 
epoch may deviate from given above. In Eq. 02]), e is a dimensionless parameter 
which determines the cutoff below which the velocity-dependence disappears. Since we 
are interested in the case that e <^ 1, we recover (av) ~ {crv)o in the hmit v — )• vq. 

The power law index n and the cutoff parameter e depend on models. The Sommerfeld 
enhancement predicts n = 1 and e ~ m^/m, while it also predicts n = 2 in the zero-energy 
resonance region. In the Breit-Wigner enhancement, the DM annihilation cross section 
reduces to the form 03]) with n = 4 and e = [6^ + 7^]/fo in the limit ^ max[(5, 7]. 
In another limit v"^ <C max[(5, 7], it is of the form with n = 2 and e = [5^ + 7^]/(25fQ), 
after rescaling {av)Q — )■ {av)ov'^/26. When the annihilation cross section with n = 4, 
the DM annihilation cross section is large enough to reduce the DM number density 
significantly even below the freezeout temperature. Thus, we consider the cases of n = 1 
and 2 in the following analysis, and it would give conservative bounds on the Breit-Wigner 
enhancement. 

With the annihilation cross section being given, the annihilation term in the Boltz- 
mann equation, which governs the evolution of the number density of DM ubm, is given 
by 



dn 



DM 



dt 



-nl^{av). (4) 



In deriving constraints from BBN and CMB, spectra of injected energy per unit time are 
needed for all the daughter particles. For the particle species i, such a quantity is given 
by 



dt 



1 2 / x^Ni 

2^dmW^, (5) 



where dNi/dE is the energy spectrum of i from the pair annihilation of DM. The energy 
spectra of decay products depend on the property of DM; for a given decay process, we 
calculate dNi/dE by using PYTHIA package |28j . 

In order for a qualitative understanding of the effects of velocity-dependent cross 
section, it is instructive to consider the total energy injection Ap in typical cosmic time, 
which is ~ H~^i with H being the expansion rate of the universe. For this purpose, let 
us define 

Ap ^ l E^;^n^y[{av)H-^ 

s ~ 2 s ' ^ ^ 
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where i^vis is the total release of visible energy in one pair-annihilation process of DM, 
and s is the entropy density. Because we consider the case that the cosmic expansion is 
(almost) unaffected by the DM annihilation, the quantity Ap/s is approximately propor- 
tional to the amount of injected energy in a comoving volume per Hubble time. Numeri- 
cally, we obtain 



^ 2.9 X 10- GeV ( ^) ( ^) (l^) | „„„ 



1 keV / \ 2m J \ m y \ 1 pb 



-, Re 
e 



(7) 



where T is the cosmic temperature. In addition, for the convenience of the following 
discussion, we have introduced the enhancement factor 



We have set the present DM energy density to be consistent with the WMAP observation, 
n^h^ ~ 0.11 [29]|1| In the case of velocity-independent cross section (where (av) = {av)o), 
it is evident that the energy injection per comoving volume decreases as T decreases. This 
is natural since the DM number density decreases as the Universe expands. In the case 
of the velocity-dependent cross section, however, some non-trivial features appear. First 
note that the velocity is estimated as 



- = J^ocTV^ for T>Tkd, 

Vq \ m 



V /25Tkd T 

(X I for I < ikd, 



(9) 



Vq \ m Tkd 

where Tkd denotes the temperature at the kinetic decoupling^ Below this temperature, 
a DM particle cannot maintain kinetic equilibrium with thermal plasma, and hence it 
propagates freely and loses its momentum only adiabatically by the Hubble expansion. 
Notice that typical kinetic decoupling temperature for WIMP DM is much smaller than 
the freezeout temperature (we may say Tkd ~ keV - MeV for WIMP DM candidates) 
[32l l33l [341 [35] . Thus it is found that, for T < Tkd and n > 1, the energy injection 
([7j) is constant, or increases as T decreases as long as the velocity dependence is not 
saturated. In Fig. [T] we plot Ap/s as a function of time for n = 1 (top) and n = 2 
(bottom). Red solid lines correspond to Tkd = 1 MeV, and green dashed lines correspond 
to Tkd = 1 keV, for e = 10~^ - 10~^ from bottom to top. Tkd = 1 keV approximately 
corresponds to a lower bound on the temperature of kinetic decoupling in order not to 
suppress the density fluctuation for a formation of the Lyman-a clouds [HE]- We have 



^ In the standard thermal rehc scenario, {crv)o is fixed once we fix the DM abundance, though it 
deviates from the canonical value ('^ 3 x 10~^^ cm'^/s) due to the velocity dependence around the DM 
freezeout epoch, especially in the case of n = 2 [30l[31]. In the following we derive upper bound on {av)o 
with fixed DM abundance. 

^Precisely speaking, there appears a dependence on the relativistic degrees of freedom g*^ for T < Tkd, 
as V (x (7*s(r)^/'^T. We have taken into account this correction. 
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taken {av)o = 3 x 10~^® cm^/s, m = 1 TeV and i?vis = 2m. Then, even at T ~ 0.1 MeV, 
we get an enhancement factor of the order of ~ lO'^ (~ 10^) with n = 1 [n = 2). 
This imphes that constraints become stronger than the case of usual DM without velocity 
dependence. 

In the following we perform detailed calculations of the effects on BBN and CMB, and 
derive constraints on the cross section with various choice of e and T^d- 

3 Constraints from BBN and CMB 

3.1 Constraints from BBN 
3.1.1 Basic picture 

It has been known that injection of high-energy particles which are emitted through 
the annihilation of long-lived massive particles during/ after the big-bang nucleosynthesis 
epoch (at a cosmic time t = 10~^ - 10^^ sec) significantly changes the light element 
abundances [18], [191 [201 [21], [221 [37] . However, the effect of the injection highly depends on 
what particles are injected. We discuss two possibilities: (i) injection of electromagnetic 
particles and (ii) injection of hadronic particles in this section. 

The injection of high-energy electromagnetic particles such as photon and electron 
induces the electromagnetic cascade, which produces a lot of energetic photons. Those 
photons destroy the background ^He and produce lighter elements such as deuterium 
(D), tritium (T), '^He, and heavier elements such as ^Li nonthermally at t > 10® sec. In 
particular there is a striking feature that the ^He to D ratio (^He/D) tends to increase. 
By comparing to the observed value of ^He/D, this gives us the most stringent bound 
on the annihilation cross section in case of the injection of electromagnetic particles |22j . 
This reaction occurs at the cosmic temperature of T ~ 10~^ MeV. It is notable that this 
constraint from BBN [38] is stronger than that on the fi- or y-distortion from the Planck 
distribution of CMB 

On the other hand, the injection scenario of high-energy hadrons such as pion, pro- 
ton {p), neutron (n) and their antiparticles might be more complicated, but has been 
understood in detail [381 HO]- The emitted high-energy neutron and proton destroy the 
background ^He and produce D, T, ^He or ®Li. The charged pions, nn and pp pairs in- 
duce an extra-ordinal interconversion between the background proton and neutron, which 
makes the neutron to proton ratio (n/p) increase. Then this mechanism produces more 
^He. In terms of the annihilating dark matter, the overproduction of D or the increase of 
^He to deuterium ratio (^He/D) gives us the most stringent constraint on the annihilation 
cross section [211 122] . 

In the following, we perform a detailed calculation of the light-element abundances; 
to take account of the injection of hadronic and electromagnetic particles, we follow the 
procedure given in [38]. Then, comparing the theoretical prediction with the updated 
observational constraints, we derive precise upper bounds on the annihilation cross section 
as a function of the DM mass. 
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Figure 1: Energy injection from DM annihilation per entropy density per Hubble time 
as a function of time for n = 1 (top panel) and n = 2 (bottom panel). Red solid lines 
correspond to T^d = 1 MeV and green dashed lines correspond to T^d = 1 keV, for 
e = 10~^, 10~^, 10"'', 10"^ from bottom to top. We have taken (cri;)o = 3 x 10"^^ cm^/s 
and 171 = 1 TeV. 
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3.1.2 Observational light element abundances 

Next we discuss observational limits on D/H and '^He/D which are adopted in this study. 
In the previous work [22], it was shown that these elements give us more stringent con- 
straints than the others. The recent observation of the metal-poor QSO absorption line 
system QSO Q0913+072, together with the six previous measurements, leads to value of 
the primordial deuterium abundance with a sizable dispersion 

(nD/nH)p = (2.82 ± 0.20) x 10"^ (10) 

Compared with data adopted in the previous analyses in Refs. f2T{ [22], the error of 
(nD/nH)p has been reduced by about 20 %. 

We adopt an upper limit on nsHe/'^D which is recently observed in protosolar clouds [12], 

(mHe/^D)p < 0.83 + 0.27. (11) 

This value was also used in Ref. [22j . 

3.1.3 Constraints on electromagnetic particle injection 

Here we discuss the case of an electromagnetic annihilation modes into electron and/or 
photon. It is notable that the total amount of energies into electromagnetic modes ap- 
proximately determines the bound, independently of the detail of each mode. In Fig. [2] we 
plot the upper bounds on the annihilation cross section obtained from the observational 
limit on ^He/D, with n = 1 (top) and n = 2 (bottom) for various values of e = 10^^" 
- 10~^. Here the kinetic decoupling temperature is set to be 1 MeV. The dashed line 
denotes the canonical annihilation cross section (= 3 x 10~^^cm^/sec ). In the top panel, 
we see that the bounds highly depend on the cutoff parameter e when e > 10~^. This 
behavior can be understood from the fact that the production of ^He becomes most effi- 
cient when T ~ 10~^ MeV; at such a temperature, the enhancement factor is estimated as 
~ 5 X 10-^(Tkd/MeV)"^/^ (m/TeV)"^/^ (T/lO-^MeV), which becomes smaller than 
~ 10~^ with the present choice of parameters. Then, when e > 10~^, the cross section 
is enhanced purely by the factor of e~^. To allow the canonical value of the annihilation 
cross section for a few TeV mass of dark matter, we need e > 10~^'^ at least. In the 
case of n = 2 which is plotted in the bottom panel of Fig. [21 R^^ is much smaller than 
e everywhere in this parameter space. Therefore determines the enhancement of the 
annihilation cross section, and there exists a simple scaling law for the line of the limits, 
which means that the upper bound is proportional to e. 

This feature is slightly different in case of Tkd = 1 keV. Because the inverse of the 
enhancement factor with = 1 is the order of ~ 1 x 10~^, any constraints with 
e < 1 X 10~^ is insensitive to e. In case of = 2, the constraint is same as the bottom 
panel of Fig. [2] because of the same reason. 

3.1.4 Constraints on hadron injection 

When we consider the injection of hadronic particles, the limit is completely different 
from that of the electromagnetic particles. The constraint on the overproduction of the 
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Figure 2: Upper bound on the annihilation cross section obtained from the observational 
■^He/D limit with n = 1 (top) and n = 2 (bottom) for various values of e = 10^^*^-10^^. 
Here DM is assumed to annihilate purely radiatively into electron and/or photon. The 
kinetic decoupling temperature is set to be 1 MeV. The dashed line denotes the canonical 
annihilation cross section (= 3 x 10~^^cm^sec~^ ). 
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Figure 3: Same as Fig. [21 but for the kinetic decoupling temperature set to be 1 keV. 
The case of n = 2 is completely same as the bottom panel of Fig. [2l 
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deuterium due to the ^He destruction often gives the most stringent constraint [22j. To 
study the hadronic injection, hereafter, we assume DM annihilates into a ly-boson pair as 
a typical hadronic DM annihilation channel; in such gnificant amount of hadrons 

are produced by the subsequent decay of the W bosons produced by the DM annihilation. 
Constraints do not change much for other cases, such as DM annihilation into bb [22]. In 
Fig. Hlwe plot the upper bound on the annihilation cross section obtained from the obser- 
vational limit on D/H with n = 1 (top) and n = 2 (bottom). The kinetic decoupling tem- 
perature is set to be 1 MeV. First let us consider the case of n = 1. Because the hadrodis- 
sociation processes become most effective at T ~ 10~^ MeV, for which the enhancement 
factor is estimated as R-^ ~ 5 x lO'^ (Tkd/MeV)"^/^ (m/TeV)"^/^ (T/lQ-^MeV), the 
constraint is determined only by the value of e if e > 10~^. To agree with the canonical 
annihilation cross section, we have to assume e > 10^^ for a few TeV mass of dark matter. 

If the kinetic decoupling occurs at around 1 keV, the enhancement factor behaves 
differently from the case of Tkd = 1 MeV because the hadrodissociation occurs before the 
time of the kinetic decoupling. As is shown in Fig. |5l then the enhancement factor is 
estimated to be R'^ ~ 5 x 10"^ (T/lQ-^MeV)^/^ (m/TeV)"^/^ at T = lO'^ MeV, from 
which we easily find that the constraint is independent of e for e < 10^^ 

When we consider n = 2, the cutoff parameter determines the upper bound everywhere 
in the current parameter space for both Tkd = 1 MeV and 1 keV. The result is shown in 
the bottom panel of Fig. H] as a representative of both cases. 

Although so far we have discussed the limit from D/H, consideration of other light 
elements sometimes tighten the constraint. In particular the limit from '^He/D by the 
photodissociation of ''He could also give us stronger limits. Note that even in the annihi- 
lation into quarks and gluons, the photodissociation occurs because a sizable amount of 
the electromagnetic particles is also injected as decay products. For example, the electro- 
magnetic energy corresponds to ~ 47% of the total energy in case of the annihilation into 
a VT-boson pair |22]. In Refs. [211 l22l |37], it has been shown that the upper bound from 
^He/D due to the photodissociation accompanied with the hadronic annihilation is much 
weaker than that from D/H when the annihilation cross section does not depend on v. 
The bound on e from D/H was severer than the one from '^He/D by three or four order 
of magnitude. On the other hand in case of the w-dependent cross section, the situation 
can be altered since the enhancement factor depends on the temperature. 

Because the constraint on ^He/D is sensitive to the cosmic history at T ~ 10"'^ MeV, 
the enhancement factor is the order of ~ 10'' for Tkd = 1 MeV with n = 1. For a small 
cutoff parameter e < 10~^, then the constraint from ^He/D can become stronger than 
that from D/H for m > 1 TeV. This feature is shown in the top panel of Fig. |6l In this 
figure we plot the constraint only from ^He/D, ignoring the D/H constraint. Notice that 
large amounts of D are produced in most parameter space as is seen from Fig. |H Since 
hadro/photo-dissociations of ^He also create ^He, both D and ^He are produced from 
the standard BBN and hadro/photo-dissociation processes. For sufficiently large {(7v)o, 
both D and ^He are produced by the hadrodissociation of '^He at T ~ 10~^MeV, and the 
constraint comes from the additional photodissociation effects at around T ~ 10~^MeV. 
These mixed processes complicate the behavior of the lines in Fig. El There is no simple 
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Figure 4: Upper bound on the annihilation cross section obtained from the observational 
D/H limit with n = 1 (top) and n = 2 (bottom) for various values of e = 10~^° - 
10^"^. Here DM is assumed to annihilate into a ly-boson pair. The kinetic decoupling 
temperature is set to be 1 MeV. The dashed line denotes the canonical annihilation cross 
section (= 3 x 10~^^cm^/sec ) which gives the right amount of the dark-matter relic 
density. 
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Figure 5: Same as Fig. HI but for the kinetic decoupling temperature set to be 1 keV. The 
case of n = 2 is same as the bottom plot of Fig. HJ 



scaling law among lines with respect to the line of f -independent constraint. On the other 
hand, if we take n = 2 , the bound from ^He/D is always weaker than that of D/H. This 
is clearly shown in the bottom panel of Fig. [61 

We also consider the case of Tkd = 1 keV. Results are shown in Fig. [71 As in the 
previous case, the constraint is mostly from the abundance of D. In some parameter 
region, however, ^He/D gives the most stringent constraint. This fact is seen in the case 
of = 1 for m > 1 TeV and e < 10~^. In addition, for n = 2 (the bottom panel of 
Fig. [7]), a simple scaling law (i.e., the proportionality of the upper bound on {av)o to e) 
breaks down once e becomes smaller than 10"*^; for such a small value of e, the constraint 
becomes significantly stringent. This feature can be understood analytically because, for 
e < 10~^, Ap/s starts to increase as a function of t after t = 10^ sec (T = 1 keV). Such a 
behavior is clearly seen as the dashed line in the bottom panel of Fig. [H 

Before closing this subsection, we comment on the constraints from the Li abundances. 
Photo/hadro-dissociation processes also modify abundances of ^Li and ^Li. However, 
constraints from observations of ''Li/H and ^Li/H are weaker than those from D/H and/or 
^He/D for annihilating DM [22]. 
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Figure 6: Same as Fig. 121 but for DM annihilating into a VT-boson pair. 
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Figure 7: Same as Fig. |3l but for DM annihilating into a Vl^-boson pair with n = 1 (top) 
and n = 2 (bottom) . 
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3.2 Constraints from CMB 



Energy injection around the recombination epoch affects the CMB anisotropy ^51 
231 [25|jl This is because injected energy can ionize neutral hydrogens and modify the 
standard recombination history of the Universe. The effect is characterized by the quantity 
z', 2;), which represents the fraction of injected electron (photon) energy E for 
i = 6(7) at the redshift z' used for ionization of the hydrogen atom at the redshift between 
z and z + dz. The evolution equation of the ionization fraction of the hydrogen atom, Xe, 
includes the following additional term, 

dz' nlyi{z'){(Tv) m dx^^lim, z' , z)_ ^^^^ 







dz 


-J 

DM J 



H{z'){l + z') nniz') E^^ dz 



where E-^j = 13.6 eV is the Rydberg energy, nn is the number density of the hydrogen 
atom and 



dxill{m,z',z) _ f 



dz J m 



'dNP d^^{E, z\ z) ^ 1 rfiVf ) dxL"^(E, z\ z) 



dE dz 2 dE dz 



(13) 



Here F denotes the final state of the DM annihilation, e.g., F = e^e , W^W , etc., 
and dNe^^ /dE denotes the energy spectrum of the electron (photon) generated from the 
cascade decay of F. In the case of F = e"'"e~, we have dNe^^ /dE = 5{E — m). For general 
final state F, it is evaluated by the PYTHIA code [2B]- We follow the methods described 
in Refs. [251 SS] to compute dxion{E,z',z)/dz. This term is included in the RECFAST 
code [17] implemented in the CAMB code [3H] for calculating the CMB anisotropy. 

Additional energy injections from DM annihilation around the recombination epoch 
cause ionization of neutral hydrogen atoms. Thus the effect is to slow down the recom- 
bination of the Universe. As a result, anisotropics in CMB are dumped at small scales 
due to the increase in thickness of the last scattering surface. Fig. |8] shows the TT power 
spectrum of the CMB temperature anisotropy, with/without DM annihilation effect. The 
solid line corresponds to the best-fit ACDM model without DM annihilation, and dotted 
line to DM annihilation cross section into e~^e~ with (av) = 10~^^cm^/sec and dashed 
line to DM annihilation cross section with (av) = 10~^^cm^/sec, while all cosmological 
parameters are fixed. Here we have taken m = ITeV with velocity-independent annihi- 
lation cross section. It is seen that DM annihilation effects suppress the TT spectrum, 
reflecting the increase in thickness of the last scattering surface. 

It is not hard to imagine that this effect has a degeneracy with other cosmological 
parameters. In particular, the increase of the reionization optical depth causes similar 
effects. In order to derive conservative bounds on the DM annihilation cross section, we 
must take into account degeneracies between DM annihilation effect and other cosmolog- 
ical parameters. We have derived 2a constraints using a profile likelihood function where 
the other cosmological parameters including the six standard ones (wb, Uc, ^a, f^-s, in 



DM annihilation also induces CMB spectral distortion, which is constrained from COBE FIRAS 
measurement ^43j . This constraint is weaker than that from anisotropy measurements [351 Slj • 
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Figure 8: Power spectrum of the CMB anisotropy with no DM annihilation effect (sohd), 
with {av) = 10~^^cm^/sec (dotted) and {av) = 10~^^cm^/sec (dashed) for m = ITeV and 
assuming DM annihilation into e~^e~ with velocity-independent annihilation cross section. 
Also shown are data points from WMAP, QUaD, ACBAR and CBI. 



the notation of Ref. |19]) and the amplitude of the Sunyaev-Zel'dovich effect are marginal- 
ized so that the original likelihood function is maximized for given DM annihilation cross 
section and mass. The likelihood surface is scanned by using the CosmoMC code [50]; in 
our analysis, we have modified the CosmoMC code to take account of the above mentioned 
effect of energy injection. The used datasets include WMAP |19], ACBAR |5l], CBI [52] 
and QUaD [53]. As opposed to BBN constraints, CMB constraint depends on the injected 
radiative energy, hence purely leptonic annihilation is more strongly constrained than the 
hadronic one. 

The result is presented in Fig. [9] where we plot upper bounds on the annihilation cross 
section obtained from CMB anisotropy data as a function of DM mass for e = 10~^ - 
10~^. DM is assumed to annihilate into e~^e~ pair in the top panel and W^W~ in the 
bottom panel. Here we have taken n = 1 and T^d = 1 MeV. We have checked that the 
results do not change for n = 2 and/or T^d = 1 keV. This is because the CMB constraint 
is sensitive to the annihilation rate at around the recombination epoch, T < 1 eV, and 
hence the annihilation cross section is already saturated for most interesting range of e for 
both n = 1 and n = 2. Comparing them with BBN constraints, it is found that the CMB 
constraint is severer for the leptonic annihilation case independently of the parameters. 

In the case of hadronic annihilation with m < a few TeV, the situation is not so simple. 
For n = 1 and Tkd = IMeV, CMB gives weaker constraint than BBN for e > 10"^, as seen 
from Fig.m but becomes tighter for e < 10~^. The situation is similar for Tkd = IkeV. This 
is because the BBN constraint from the observation of D /H is sensitive to the annihilation 
at T ~ 10~^MeV, and the annihilation cross section do not saturate at that epoch for 
small e for n = 1. On the other hand, for n = 2, BBN gives tighter constraint than the 
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CMB for parameter ranges shown in the figures. 

Therefore, CMB takes complementary role to BEN in constraining the DM annihila- 
tion with velocity-dependent annihilation cross section. 

4 Conclusions 

In this paper we have investigated effects of DM annihilation on BEN and CME, and 
derived constraints on the DM annihilation rate, particularly focusing on the case where 
the annihilation cross section has a velocity- dependent structure. This is partly motivated 
by the observations of cosmic-ray positron/electron excesses and their explanations by 
the DM annihilation contribution. We phenomenologically parametrized the velocity- 
dependence of the annihilation cross section and the critical velocity at which such an 
enhancement saturates, and derived general constraints on them. Our constraints are 
applicable to known velocity-dependent DM annihilation models, such as the Sommerfeld 
and Ereit-Wigner enhancement scenarios. 

These results have been plotted in Figs. |l]-|n]by changing the parameters and obser- 
vation. Therefore readers can read off the allowed parameter regions from those figures 
in accordance with the intended use. 

Some comments are in order. If the DM annihilation is helicity-suppressed, the p-wave 
process may be the dominant mode, as is often the case with Majorana fermion DM. In 
this case we obtain n = —2 : (av) oc v"^. Thus the annihilation cross section becomes 
smaller as the temperature decreases, until the S'-wave process becomes efficient. For 
negative n, the EEN/CME constraints are weaker than the velocity-independent case. 

In the Sommerfeld enhancement scenario, it was pointed out that the DM-DM scat- 
tering mediated by light particle exchanges causes observationally relevant effect [51], and 
this also gives significant constraint . 
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